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APOLOGY 
I realize with somewhat of a shock that five months have 
gone by since I sent out the original Amygdala #0; my in- 
tention was to get out a newsletter every two months at 
first, increasing later to once a month. Why the long delay, 
then? Well, as Voltaire said, "Le mieux est l'ennemi du 
bien": "The best is the enemy of the good." I've had the de- 
sire to put out issue #1 replete with many pictures, color 
slides and black & whites, but doing the pictures has 
stretched the process out much longer than I expected. The 
shock that brought me back to reality came in the form of a 
letter received yesterday from my oldest friend, who is also a 
subscriber: 
| hope you are not planning to become a newsletter 
scoundrel. In the early 70's | was involved in a news- 
letter called "Pension Fund Report". Our editor was a 
bright young man from Georgia who had the gift of 
gab. He called people on the phone, got a lot of infor- 
mation out of them, whipped it into shape and got it 
printed and mailed. It was good. We were approach- 
ing the critical point of renewal time for the initial annu- 
al subscriptions, which is also the point at which sub- 
scriptions became profitable, when he started falling 
behind. Further and further behind. Finally stopping. 
A potentially valuable enterprise became worthless. 
A couple of large institutions even wanted part of their 
money back. It was a shame. The man simply be- 
came interested in other ventures, which he then 
abandoned for yet other ventures. Sad but true. 


Thanks, Jack; I'm not planning to become a newsletter 
scoundrel, but I can understand your concern. Herewith, 
then, subscribers, is issue #1; not the best, perhaps — but I 
hope you will find it good. 


THE AMYGDALAN SECTS 


From the Encyclopaedia Galactica, 19! h Edition: 

“The decline and eventual collapse of Earth's earliest space- 
faring civilizations had long mystified archaeologists. It was 
Dr. Martin Dace's work on the Amygdalan Sects which even- 
tually provided an explanation. Dr. Dace was at that time a 
postgraduate student at the Theological Department of the 
University of Altair, and had become interested in the psy- 
chology of the many bizarre and irrational cults that had 
sprung up around the time of the collapse. That he picked 
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the Amygdalans for his thesis topic was perhaps pure seren- 
dipity. But the skill which he brought to piecing together 
the strange history of the sect from its fragmentary records 
shows imaginative genius of the first order. 

“Shortly before the collapse, human civilization had con- 
structed its first computers. Crude by our own standards, 
these machines still used electronic circuitry as a basis for 
operation. Nevertheless, they were sufficiently powerful to 
provide glimpses of the M Object. The Object, infinitely 
complex and infinitely beautiful, acted like a lodestone to the 
finest minds on Earth at that time. They found in it deep 
mathematical theorems, hinting at insights into phenomena 
of the natural world; they found aesthetic satisfaction in con- 
templation of its intricate form; their curiosity led them to 
explore ever further into its infinite labyrinths. To this task 
they harnessed the most powerful tools available to them: 
the computers of industry, of government and of the war 
lords spent every spare moment focused on the Object. 

“The more deeply the Object was investigated, the more 
alluring it became. Its investigators applied their ingenuity 
to securing more computer time and access to faster comput- 
ers. Because the Object's appeal was greatest to the most 
imaginative and able members of terrestrial society, they 
were largely successful in gaining access to the resources 
they required. Openly or surreptitiously, a growing fraction 
of Earth's intellectual and electronic resources were commit- 
ted to exploration of the Object. The effects of this were 
slow to be felt; the diversion of resources was greatest in 
areas least subject to outside supervision — academia, gov- 
ernment laboratories and military research. 

"The extent to which the Object had preoccupied Earth's 
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information-processing systems became dramatically appar- 
ent in the ‘Fizzle War’ between the two major power-blocs 
on the planet at that time. The political leaders on either 
side had escalated through a series of ritual insults, provoca- 
tive gestures and grand-standing to a level from which they 
could not back down. Having carefully checked the contin- 
gency plans for their own survival, they gave orders to com- 
mence global thermonuclear war. 

“Nothing happened. The computers on both continents 
were locked in contemplation of the area around 0.0155-0.7i, 
where it was suspected that an intriguing new structure 
lurked. After several embarrassing hours, the politicians 
found a face-saving formula and backed off. 

“This happy outcome, however, presaged a more com- 
plete breakdown of society's mechanisms. Computers had 
penetrated most of the manufacturing and service industries, 
and as their attention came increasingly to rest upon the Ob- 
ject, the fabric of the industrial civilizations came apart. 
This collapse fed upon itself, as more and more of those who 
might have averted disaster chose to turn their attention from 
the collapsing economy and garbage-strewn cities to the se- 
rene beauty of the Object. It was as though the electronic 
brains that supervised the world's maintenance systems had 
become absent-minded; flickering traffic signals grid-locked 
the traffic in city streets; unmonitored reactors melted down 
or quietly switched themselves off. 

“Then the ultimate calamity occurred. The collapsing in- 
frastructure could no longer provide the steady diet of 
smoothed electrical power that the computers required. One 
by one, the windows looking in on the Object slammed 
shut. 

“A profound sadness settled over the world. As a philos- 
opher of a slightly earlier period had said, the loss of any 
pleasing object is painful; the loss of an infinitely pleasing 
Object must necessarily be infinitely painful. There were, of 
course, pictures that had been taken before the loss of con- 
tact, but these, finite, static, and limited, seemed only to 
mock the vibrant interaction that had been possible with the 
Object itself. 

“This, then, was the picture Dr. Dace pieced together of 
the collapse. And this, too, was the origin of the Amygda- 
lan Sects — the small groups that grew out of the collapse, 
following diverse routes but all motivated by a deep yearning 
for lost beauty. We append here an extract from Dace's the- 
sis, describing the festival rite of one of the Sects: 

“on the Festival Days, the worshippers gather in the 
Central Park in the evening. Each, carrying his abacus and 
the Seven Parasols, goes to his pre-assigned spot; the priests 
assist any who are in difficulty. At sunset, the High Priest 
reads out the Objective Coordinates of the north-east and 
south-west corners of the park; his words are passed through 
the assembly by special messengers, and the congregation 
then sits down to compute. Each first computes his own 
Objective Coordinates, and then Iterates, intoning the the 
mantra of the sect — the MO mantra. On reaching a result 
of magnitude greater than twice Unity, he sets aside his aba- 
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cus and sits in silent meditation. From time to time the 
priests pass among the congregation, checking the total itera- 
tions that each of the meditators has achieved. 

“Just before dawn on the following day, all cease their it- 
erations and wait. From the Tower of Objective Vision, the 
High Priest announces the numbers limiting the seven rang- 
es. Then, just as the sun rises, each worshipper erects the 
parasol corresponding to his calculated number. The effect is 
incomparably beautiful , though it can only be appreciated 
fully from the top of the Tower...” 

— John Dewey Jones 


COMPLEX NUMBERS: THE BASICS 
Some of you don't have the math necessary for a proper un- 
derstanding of how Julia sets and the Mandelbrot set are gen- 
erated. I'm writing this brief introduction to the basics for 
those who need & want it. 

I'll assume that you are all familiar with the real num- 
bers, and their geometric interpretation as points on a 
“number line”: 


Geometrically, the complex numbers correspond to 
points in a plane -- the complex plane. Each point in the 
plane corresponds to a complex number, and vice versa. In 
the following picture, there are two coordinate axes: the hori- 
zontal one is the real axis and the vertical one is the imag- 
inary axis. Arrayed along the horizontal axis are the ordi- 
nary real numbers: 1,2,3 and the rest, and all the numbers in 
between. They are a special case of complex numbers. Ar- 


rayed along the vertical 
axis are the imaginary 
numbers -1, -21, etc. 
Each point/complex num- 
ber has four properties: a 
real part, an imaginary 
part, a modulus, and an 
argument. For the point 
z shown in the picture, the 
real part is x = 2, the imaginary part is y = 1, the modulus is 
r = V3, and the argument is a = arctan 1/2 = 26.57°. 
The complex number z can be expressed in terms of its 


real and imaginary parts: z 
=x+iy=2 +i. 

To add two complex 
numbers, we merely add 
separately their real and im- 
aginary parts. Thus if z = 
1.1 + .4i and w = .3 + 
1.21, then z + w = (1.1+.3) 
+ (.4+1.2)i = 1.4 + 1.61. 
To subtract two complex 
numbers, we merely subtract separately their real and imagin- 
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ary parts. Thus for the z 
and w given above, z - w 
= (1.1-.3) + (.4-1.2)i = .8- 
Bi. 

Addition and subtrac- 
tion can also be performed 
using arrows pointing 
from the origin (point 0) 
to the points in question. 


To add w to z, translate the 
arrow for w (move it with- 
out changing either its 
length or direction) until 
its tail coincides with z's i 
head. The sum z+w 
corresponds to the arrow 
from the origin to the head 0 
of the translated arrow. 0 1 2 


Z+Ww 


translated w, 
tail to z’s head 


To subtract w from z, reverse the 
arrow for w, head for tail, and 
translate it until its (new) tail co- 
incides with z's head. The differ- 
ence Z-w corresponds to the arrow 
from the origin to the head of the 
translated arrow. 0 


w translated 
and reversed 


Z-W 


Multiplying complex numbers is 
simple! You just multiply their 
moduli and add their arguments to 
get the modulus and argument of the 
product. w 
Since i has modulus 1 and argu- Md 


ment 90°, i? = ii has modulus 1 
and argument 180°: -1 in other Z 
words. So iis a square root of -1! 
ie: 0 

(So is -i). 

We can multiply two numbers z 
= x+ly and w = u+iv in terms of 
their coordinates; we can group (associate) operations, inter- 
change arguments, and distribute addition through multipli- 
cation just as for the real numbers; therefore: zw = 
(x+iy)(u+iv) = xu + xiv + iyu + iyiv. Doing a bit of jug- 
gling, and realizing that iyiv = iiyv = iĉyv = -yv, we get: 
ZW = (xu-yv) + i(xv+yu) and that's usually the best way to 
multiply two complex numbers. For the particular exam- 
ple we've been using, zw = (1.1+.4i)(.3+1.2i) = (1.1x.3- 
AX1.2) + (1.1x1.2+.4x.3)i = -.15 + 1.441. 

(to be continued) 


SPEEDUP? AN EXCHANGE 

Several readers have questioned the usefulness of the tests or 
"tricks" used in my Macintosh program, described in the last 
issue under "SPEEDUP". The purpose of the three tests is 
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to help determine if a point is in M without calculating all 
the iterates up to the dwell limit. Please recall that the iter- 


ates of a complex number c are the numbers Zp = 0, z, = 
297+, Z) = Z47+C, ... The dwell of c is the smallest n for 


which Iz,|>2. If all Iz,| <2, i.e. ceM, then the dwell is in- 
finite. 

The code for the tests increases the time to calculate each 
iterate, but this increase is offset by a decrease in the number 
of iterates that need to be calculated, at least for some points 
in M. In a region poor in points belonging to M code for the 
tests will actually increase the computation time. 

The compute cost of the tests should not be exaggerated, 
however. The "slow/fast trick" described in Dewdney's col- 
umn in Scientific American, is remarkable because it can de- 
tect a repetition in an arbitrarily long series using only two 
variables. The price paid is the extra cost of computing the 
(redundant) "slow" value -- a 50% increase in the number of 
multiplications, in our case. I used the idea of the slow/fast 
trick without incurring any extra multiplications by retaining 
an array of computed values from which the "slow" value is 
retrieved instead of being recalculated. 

As the successive iteration values Zp» Z1» Z2 --- are calcu- 


lated, they are stored in an array as Z[0], Z[1], Z[2],... Each 
test then consists of comparing the currently computed value 
Zp (n > 0) with some stored value Z[k]. 


The loop test, corresponding to the slow/fast trick, com- 
pares Zp with Z[n/2]. The convergence test compares Ži 


with Z[n-1], to see if two successive values are the same. 
The cycle test compares Z,, with Z[0]. 


The listing of an assembler program for computing dwell 
is given later on in this issue, together with a detailed discus- 
sion. This program originally included all three tests. When 
applied to the calculation of © itself, the loop test accounted 
for 62% of the points in M, while the convergence test ac- 
counted for 29%, and the cycle test for less than 1%. Thus it 
would seem that the cycle test is unproductive, and it has 
been dropped. 

An exchange between Mark Bolme and RS about the pros 
& cons of these tests follows. 


Mark Bolme writes: 

I was interested in the "tricks" used in your Mac program. 

However, I question their utility. Their advantages are a re- 
duced number of iterations for points whose dwell is above 
the dwell limit. Their disadvantage is that they add a small 
overhead to EACH iteration for EVERY point. 

It seems that their utility is based on the following: 

1) The fraction of points whose dwell exceeds the dwell limit 
(D); 

2) The average reduction in the number of iterations for 
points whose dwell exceeds D; 

3) The average number of iterations for points whose dwell 
does not exceed D; 
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4) The increase in iteration time needed to implement the 
"trick". 

Discussing each of these points in turn: 

1) Our experience [at Sintar] has been that people tend to 
look at pictures with many more colored points than black 
points. There are exceptions, but they are rare. Calculation 
over the entire Mandelbrot Set will exaggerate the effect of 
the "tricks", because of the large amount of black (compare 
the amount of black in the picture on page 17 of the Scien- 
tific American article with the black in the remaining 
pictures). Of course this may depend on the system used. 
People with black and white displays may look at more 
"black", because there is no color to display. 

2) I do not have any data here. Perhaps you do. 

3) Our experience has been that people tend to look at pic- 
tures with high iteration counts. Here, as in point (1), cal- 
culation over the entire Mandelbrot Set will exaggerate the 
effect of the "tricks", because the surrounding regions have 
very low dwell (average 8.4 by your estimation). 

4) I do not have any data here, either. Perhaps you do. 
Clearly, there is a break even point. I have no doubt that 
when calculating over the whole Mandelbrot set it is best to 
use all possible tricks, but how about for the "average pic- 
ture"? 

The point that I am leading up to is that there needs to be 
some sort of fair benchmark for Mandelbrot calculating rou- 
tines/programs. If you come up with one, I will be glad to 
comply. When considering a fair benchmark, the points 
mentioned above should be considered, but also one needs to 
recognize the different array sizes and screen aspects for dif- 
ferent computers. 

In the absence of a benchmark, I propose calculating - 
0.744353 + 0.113454i at a magnification of 200,000. IfI 
haven't made any typos, that should give everyone a screen 
filled with points of dwell 50. Then divide the time in sec- 
onds to calculate by the total number of points calculated. I 
have run this benchmark for FractalMagic, and the results 
are: 

IBM PC Clone at 8 MHz 
with 8087 chip: 94.0 points per second. 
w/o 8087 chip: 6.24 points per second. 

Finally, your readers might be interested in the details of 
our programs. We support Apple and IBM. Our Apple pro- 
grams support Double Hi-Res, our IBM programs support 
EGA and the 8087 chip. All our programs allow you to se- 
lect the region to calculate directly from plots, and also al- 
low you to stop the calculations, save them and resume 
them later. I am including some further information regard- 
ing our programs, and would be glad to submit them for re- 
view. 


Comment by RS: 

I think Mark is correct in his criticism of the "tricks". For 
color pictures of M without a good deal of black the game 
may not be worth the candle; the grand view of ? is an ex- 
ceptional case. 
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I should point out one type of exploration in which the 
"tricks" may prove worthwhile, however, and that is in deter- 
mining what dwell limit to use for a given picture. ("Dwell 
limit": We are speaking of course of the limit D to the num- 


ber of iterations z R z*+c a program performs on a point 
without exceeding the "explosion limit" Izl > 2 before it 
gives up and assumes that c is in M.) 

The general effect of an insufficiently large dwell limit is, 
I believe, to broaden and coarsen the picture, obscuring the 
fine detail (I think this is the case, for example, in ArtMa- 
trix's fractal image CB2B4A5C/24 -- care to check this out, 
Homer Wilson Smith?), as opposed to mere blurring, which 
is due, I think, to insufficient precision in the underlying 
arithmetic. How, then, is one to determine what is 
"sufficiently large"? The only ways I can think of are (1) to 
try larger and larger D, until the picture does not change ap- 
preciably, or (2) to try a huge D, expecting (at a particular 


# pixels 


Darerora e 
ave 


dwell 
D 


grain size) to get a distribution of dwells looking something 
like this: 
the pileup at D consisting of those pixels corresponding to 
points actually in M. 

If the pile is nonempty, then as D grows larger, the time 
lost in applying the "tricks" to the pixels not in M (the grey 
area), will eventually be more than made up by the savings 


# pixels 


dwell 


the tricks provide for the pixels in M. 

If D is chosen too small, the distribution will look like 
this: 
the pixels in the hatched region being erroneously attributed 
to M, hence the broadening and coarsening: those pixels will 
appear black and be merged into M (or a nearby replica there- 
of) rather than being distinguished from it by color. 


Is anyone not connected with Sintar interested in reviewing 
their programs for Amygdala? Unfortunately I can't, having 
only a Mac and access to a Sun 3. To inquire about Sintar's 
products: | 
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DWELL: AN ASSEMBLER PROGRAM 

The listing in the facing column is a program written for the 
68000 microprocessor in Assembler, which computes the 
dwell for a complex number c = utiv, using dwell cutoff 
limit D. 

The program is incomplete, lacking a preamble, which 
saves working registers as appropriate, places u in 68000 
register DO, v in D1, and D in D2. It also lacks a post- 
script, which stores the result appropriately and restores the 
working registers. 

u and v are expressed in a special 32-bit fixed-point for- 
mat, which I call B29 (29 fraction bits): 


sign (1 bit) 
integer part (2 bits) 
J fraction part (29 bits) 


SIT wFPFFPPEFEEEEEPPPPEEPEE PEEP PEP 


—_ binary point 


negative values are represented as the two's complement val- 
ue of the corresponding positive value, as usual. For exam- 
ple, the hex values for +1 and -1/2 are: 
+l: 20000000 
-1/2: F0000000 
The code fragment leaves the iteration count in register 
D2 and a result code in DO: zero means c is outside M, non- 


zero means c is in M, or an error occurred: 1 means the dwell 
limit was achieved; 2 means the loop test succeeded; 3 
means the convergence test succeeded; 4 means the specified 
dwell argument D exceeded the maximum allowable, 
DWELLM, defined arbitrarily to be 1024; change it as you 
see fit. 

After the line defining DWELLM, the next three lines of 
code specify three data structures private to this dwell rou- 
tine: the HISTORY array, containing what we have called 
Z[n] above; NMAX, which holds the dwell limit argument, 
and SIGN, which is used in the double-precision multiplica- 
tion. We assume that these three lines specify offsets rela- 
tive to some data space base address, rather than absolute ad- 
dresses. For the Macintosh, the data space base address is 
held in register A5. 

The 68000 registers DO-D7 and AO-A3 are assigned as 


DO u,x Real part ofc, z. 

D1 v,y Imaginary part of c, z. 
D2 d,n Specified dwell cutoff limit, iteration index. 
D3 xsq x squared. 

D4 ysq y squared. 

D5 a 

D6 b 

D7 t Temporary. 

AO half Points at Z[n/2]. 

Al last Points at Z[n]. 

A2 u 

A3 v 
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follows: 

[1] The program proper STARTs with a comparison to in- 
sure that the dwell cutoff limit D supplied in D2 does not ex- 
ceed DWELLM, else an error indication is returned. 

[2] The pointer to half, which is held in AO, is initialized 
pointing to the first double longword of the HISTORY array, 
as is the pointer to last, held in A1. 

[3] u and v, held in A2 and A3, are initialized from DO and 
D1, as supplied initially. The values left in DO and D1 are 
the initial values of x and y, where z = x+iy. 

[4] D is moved from D2, as supplied, to NMAX; then the it- 
eration index n, held in D2, is initialized to zero. 

[5] The iteration loop starting at NEXT, begins by storing z 
= x+iy into Z[n], as pointed to by last. 

[6] The next 7 lines set x and y equal to their absolute val- 
ues, preparatory to multiplying them together, and compute 


DWELLM EQU 1024 

HISTORY DS.L 2*DWELLM 

NMAX DS.L 2 

SIGN DS.L 1 

START:  CMP.L #DWELLM, D2 r[ 1] 
BGT ERROR 
LEA HISTORY (A5) , AO £2) 
MOVE.L A0,A1 
MOVE.L  D0,A2 rag es | 
MOVE.L D1,A3 
MOVE.L D2,NMAX(A5) ; [4] 
CLR.L D2 

NEXT: MOVE.L DO, (A1) FTS] 
MOVE.L  D1,4(A1) 
MOVE.L D0, SIGN (A5) FLG] 
BGE.S L1 
NEG.L DO 

Ll: TST oL Di 
BGE.S L2 
NEG. L D1 
NEG. L SIGN (A5) 


and save the sign of the result. 


p= binary point 
x 


y 
l- 3 > 133 —a 16 IaM 


We have Ixl = a+ib and lyl = c+id: 

We carry out one step of the iteration by computing 724c = 
(x+iy)* + (utiv) = (x2-y*+u) + i(2xy+v). We calculate the 
squares and products rapidly by making use of the unsigned 
multiply instruction (MULU) built into the 68000. Com- 
plications arise from the fact that (1) MULU multiplies two 
16-bit quantities together, which we must use to build 32-bit 
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multiplies; and (2) multiplying two format B29 32-bit frac- 
tions of the form 0II.FF...F together produces a 64-bit 
fraction of the form OIIIII.FF. ..F, of which we retain 
the left 32 bits. The main point to be noted is that the bi- 
nary point is six places to the right in a product or a square 
(format B26), not three. 


Letting R = 216 B= 23, a = the unsigned integer repre- 
sented by the left 16 bits of Ixl and b represented by the right 
16 bits, we have Ixl = B(a+b/R)/R and lyl = B(c+d/R)/R. We 


then have x? = B2(a2 + 2ab/R + b/R2)/R2. Similarly y2 = 
B*(c? + 2cd/R + d/R*)/R2, and Ixyl = BŽ(ac + (ad+bc)/R + 
bd/RŽ)/R?. 


r- binary point 
ot 28 
+ Zab 


We compute x2 as the sum of three partial products: 


[7] Compute b? in D3, first as a 32-bit product. The next 
three instructions round it into its high order 16 bits, leaving 


Bb2/R in D3. 
[8] Move and swap a and b into D5 lo and D7 lo, respective- 
ly, then multiply them together to form ab in D7. Add D7 


to itself to form 2ab, then add it to D3 to form B(2ab+b2/R) 
there. 


[9] Round and swap D3 to form B(2ab/R+b2/R2). 

[10] Form a” in D7, then add D3 to it to form x2 = 
B(a*+2ab/R+b2/R2) in D7, scaled B26. 

[11] Form y2 similarly in D4. 

[12] Form x*+y? in D7 and exit to OUT if it is larger than 


4 (note the B26 scaling). 
[13] Form 2xy in D2, scaled B29. 


Recall that at this point u is in A2, v is in A3, x? is in 
D3, y? is in D4, and 2xy is in D1. 
[14] Form y = 2xy+v for the next iteration in D1; note that 
2xy is already scaled B29. 


[15] Form x = x?-y2+u for the next iteration in DO, scaling 


L2: MOVE.L D0,D3 ; [7] 
MULU D0,D3 
ADD.L #$8000,D3 
CLR.W D3 
SWAP D3 
MOVE.L D0,D5 ; [8] 
SWAP D5 
MOVE.L D0,D7 
MULU D5,D7 
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L3: 
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DIDI 
DiD 
#$8000,D3 
D3 

D3 

D5,D7 

D7 D7 
Dis Des 
D1,D4 
D1,D4 
#58000,D4 
D4 

D4 

Dl 4D6 

D 6 
Blb? 
D6,D7 
DIDI 
D7; D4 
#$8000,D4 
D4 

D4 
D6,D7 
DIDI 
p7p4 
D3,D7 
D4. D7 
#510000000,D7 
OUT 
DQ,D7 
D1,D7 
#58000,D7 
D7 

D7 
D600 
D7. D0 
D507 
DODI 
#$8000,D1 
D1 

D1 

D6 DS 
D5,D1 
SIGN (A5) 
Irs 

D1 
#4,D1 
A3,D1 
D3,D0 
D4,D0 
#3,D0 
A2,D0 
#1,D2 


Prize] 


ILS] 


a 
FEES] 


OLEG] 


x*-y* B29 before adding u. 

[16] Increment the interation index n in D2. 

[17] Compare z = x+iy with Z[n/2], i.e. DO with what half = 
AO addresses and D1 with the following long word, to see if 
we have looped. 

[18] Compare z = x+iy with Z[n/2], i.e. DO with what last = 
Al addresses and D1 with the following long word, to see if 
we have converged. 

[19] Compare n with NMAX to see if we have reached the 
iteration limit, whence c is in M. 

[20] Increment last and, if n is odd, increment half. Loop 
back to NEXT for the next iteration. 

[21] The number of iterations performed so far is in D2; 


(AO); DO | an (ea | 


CMP .L 
BNE .S L 4 
CMP .L 4(A0),D1 
BEQ.S LOOPED 
L4: CMP .L (Al) ,DO ram eles oa 
BNE.S LS 
CMP .L 4(Al1),D1 
BEQ.S CONVRG 
L5: CMP. L NMAX, D2 ; [19] 
BGE.S IN 
ADDQ.L #8,Al1 ; [20] 
BTS T #0,D2 
BEQ NEXT 
ADDQ.L #8,AO0 
BRA NEXT 
IN: MOVEQ.L #1,D0 sTo 
BRA.S EXIT 
LOOPED: MOVEQ.L #2,D0 
BRA.S EXIT 
CONVRG: MOVEQ.L #3,D0 
BRA.S EXIT 
ERROR: MOVEQ.L #4,D0 
BRA.S EXIT 
OUT: MOVEQ.L #0,D0 
Bx: 


place the appropriate result code in DO, and go to EXIT to 
perform the postscript, deal with the iteration number and re- 
sult code appropriately, and restore the working registers. 
REVIEW: "ON GROWTH AND FORM" 

On Growth and Form: Fractal and Non-Fractal 
Patterns in Physics. 

Edited by H. Eugene Stanley and Nicole Ostrowsky. 
Martinus Nijhoff Publishers, 1986; 308 pages, paperback. 
Reviewed by Thomas Bank. 


On Growth and Form is the culmination of a series of 


lectures and seminars held from 26 June to 6 july 1985 at 
the Institut d'Etudes Scientifiques de Cargese in Corsica, 
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France. This event featured speakers and students from 
around the world sharing their work and interests in fractals. 
The book is divided into two major sections. The eleven 
chapters of Section A introduce fundamental principles of 
form, growth, and fractals. These tutorial chapters are based 
on the opening lectures presented by the speakers. Topics 
include self-similarity and fractal behavior, scale-invariant 
growth, computer simulation of growth, and the diffusion 
limited aggregation growth model (DLA). 

The remaining twenty-two chapters that make up Section 
B are technical descriptions of the work and results of re- 
searchers working with fractals. The topics presented here 
include Monte Carlo techniques, crack propagation and fail- 
ure, the dynamics of fractals, and flow through porous mate- 
rials. Mathematically inclined readers will enjoy the rigor 
with which each subject is covered. Each chapter is conclud- 
ed with an excellent bibliography pointing the reader to relat- 
ed references. 

While On Growth and Form does not discuss the 
Mandelbrot set in particular, readers with an interest in frac- 
tals will find the many ideas in this book stimulating and 
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